On the distribution of exponential functional for Levy 
processes with jumps of rational transform 
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Abstract 

We derive explicit formulas for the Mellin transform and the distribution of the exponential 
functional for Levy processes with rational Laplace exponent. This extends recent results by Cai 
and Kou [3] on the processes with hyper-exponential jumps. 
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1 Introduction 

Assume that A is a Levy process and e(g) is an independent random variable, which is exponentially 
distributed with parameter q > 0. The exponential functional of A is denned as 

e(g) 

/, = J e Xt dt. (1) 
o 

We can extend this definition to the case q = 0, if we interpret e(0) = +oo and assume that the process 
A drifts to — oo. Our main object of interest is the probability density function of I q , defined as 



p( x ) = —¥(L < x), x>0. 
ax 



Exponential functional of a Levy process is a very interesting object, which has many applications in 
such areas as self-similar Markov processes, branching processes and Mathematical Finance. See [1] for 
an overview of this topic. There exists an extensive literature covering the asymptotic behavior of p(x) 
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as x — > +00 (see [5, 6, 13, 18, 19] and the references therein) or as x — > + (see [15, 2]). At the same time, 
the distribution of the exponential functional is known explicitly only in a few special cases: when X is a 
standard Poisson process, Brownian motion with drift, a particular spectrally negative Lamperti-stable 
process (see for instance [5, 11, 16]), spectrally positive Levy process satisfying the Cramer's condition 
(see for instance [17]). Note that in all above cases we have processes with one-sided jumps, and until 
very recent time there were no known examples of processes with double-sided jumps, for which the 
distribution of the exponential functional is known explicitly. However this situation has changed in the 
last several years. 

First of all, in a recent paper [3] (which is the main inspiration for our current work), Cai and Kou 
have obtained an explicit formula for the Mellin transform of I q in the case when X has hyper-exponential 
jumps. While the authors did not provide a formula for the probability density function p(x), it does 
follow from their results rather easily with the help of the theory of Meijer's G-function. Second, in [10] 
Kuznetsov and Pardo have derived an explicit infinite series representation and complete asymptotic 
expansions for p(x) in the case of hypergeometric Levy processes. 

In this paper we will pursue the following three goals. First of all, we will generalize the results of 
Cai and Kou [3] to the class of Levy processes with jumps of rational transform, and we will also cover 
the cases when X has zero Gaussian component and/or zero linear drift. Second, we will show how to 
use the verification technique which was developed in [10]; we feel that this technique is of independent 
interest as it considerably simplifies the derivation of many results on exponential functionals. Finally, 
we will present explicit formulas and complete asymptotic expansions for p(x), which allow very simple 
numerical evaluation of this function and other related quantities, such as the price of an Asian option 
with exponentially distributed maturity. 



2 Results 

We will work with a Levy process X with linear drift and diffusion coefficients p. G R and a > 
respectively, and for which the density of the Levy measure is given by 

J raj j rhj 

tt{x) = 1 {X>0} Y Y n "''' ''''' + 1 {^<o} Y Y ^ax^e^. (2) 

j=l i=l j=l i=l 

We assume that rrij, rhj £ N, Re (p., ) > and Re(/)j) > and that pi ^ pj and pi 7^ pj for all i ^ j. It is 
clear that X has jumps of finite intensity A, where 

\ = I n(x)dx = YY - ^-pj' + Y Y - 

K j=l i=l j=l i=l 

in particular the process X can be identified as 

N t 

X t = aW t + fit + Y^ 

i=l 

where iV is a Poisson process with intensity A and £j are iid random variables, independent of N and 
with distribution P(£j e dx) = \~ l n(x)dx. 
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Formula (2) and the Levy-Khinchine formula imply that the Laplace exponent ip(z) = lnE[exp(zXi)] 
is a rational function, which has the following partial fraction decomposition 



, . . <? ■> Odjil — 1)! OLiAl — 1) , . , 

^ ) = _^ + ^ + gg_^_L + gg_A_L_ A . (3 ) 

We see that ip(z) has poles at {pj}i<j<j and {— p J } 1< < j with the corresponding multiplicities and 



j. Note that V'O 21 ) has M = Y2j=i m j P°l es (counting with multiplicity) in the half-plane Re(z) > 



and M = J2j =1 rhj poles in the half-plane Re(z) < 0. 

The poles of ip(z) are important objects for our further results, but even more important are the 
zeros of ip(z) — q, or, equivalently, the solutions to the equation if)(z) = q, where q > 0. First of all, let 
us find the total number of these solutions. Since ip{z) is a rational function of the form (3), we can 
rewrite it as 

M 

The degree of the numerator {denominator} will be denoted as Q = deg(Q) {P = deg(V)}. Note that 
(3) implies that P = M + M and 

Q = P + 2if(j>0, 
Q = P + 1 if a = and fi ^ 0, 
Q = P if a = fi = 0. 

We claim that in all three cases the equation ip(z) = q has exactly Q solutions ( = ((q) in C (counting 
multiplicity). This is obvious in the case a > or a = and /i ^ 0, as in this case the equation if}(z) = q 
is equivalent to the polynomial equation Q(z) — qV(z) = 0, which has degree Q. This statement is also 
true in the case a = ji = 0: here we have if}(z) = Q(z)/V(z) — > —X as z — > oo, therefore the leading 
term of the polynomial Q(z) — qV(z) can not be zero, thus the degree of this polynomial is also equal to 
Q (note that the degree of this polynomial is equal to Q — 1 if q — —X < 0). 

Let us denote the zeros of ip(z) — q in the half-plane Re(z) > {Re(z) < 0} as £i, (2, ■ ■ ■ , Ck 
{— Ci, — C2, • • • , — Cr}- Note that by definition we have Q = K + K. We also assume that these numbers 
are labelled in the order of increase of the real part, so that Re((y) < Re(Cj+i) and Re(Cj) < Re(Cy+i). 

Let us summarize the definition of some important quantities, which will be used frequently later: 

• M {M} is the number of poles of ip(z) in the half-plane Re(z) > {Re(z) < 0}, 

• K {K} is the number of zeros of ip(z) — q in the half-plane Re(z) > {Re(z) < 0}, 

• P — M + M {Q — K + K} is the total number of poles {zeros} of ip(z) — q. 

We collect all important information about the solutions of the equation i/)(z) = q in the next 
Proposition. 

Proposition 1. Assume that q > 0. Then 

(i) Ci and (j are real positive numbers. Moreover, for all j ^ 1, £1 < Re(£j) and £i < Re((j). 
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(%%) If a >0 then K = M + 1 and K = M + 1. 

(Hi) lfa = 0,andn>0{fi< 0} then K = M + 1 and K = M { K = M and K = M + 1} . 

(iv) lfa = /j = 0, then K = M and K = M. 

(v) There exist at most 2Q — 1 complex numbers q such that equation ip(z) = q has solutions of 
multiplicity greater than one. 

(vi) As q — > + we have the following possibilities: 

if E[Xi] > 0, then Ci(0 + ) = and Ci(0 + ) > 0, 
if E[Xi] < 0, then Ci(0+) > and Ci(0 + ) = 0, 
if E[Xi] = 0, then Ci(0 + ) = and Ci(0 + ) = 0. 

Proof. The proof of (i)-(iv) follows from Lemma 1.1 in [12] (a somewhat simpler proof of a similar result 
has appeared recently in [8]). The proof of (vi) follows easily from the fact that if)(z) ~ ELYi]z + 0(z 2 ) 
as z 0. 

Let us prove (v). Assume that equation tf)(z) = q has a solution z = z of multiplicity greater than 
one, therefore ip'(z ) = 0. Using (4) we find that ip'(z ) = implies Q'(z )V(z ) — Q{z )V(zq) = 0. The 
polynomial H(z) = Q'(z)V(z) — Q(z)V'(z) has degree at most Q + P — 1 < 2Q — 1, thus there exist at 
most 2Q — 1 distinct points Z\. for which ip'(zk) = 0, which implies that there exist at most 2Q — 1 values 
q, given by qu = ip{zk), for which the equation ip{z) = q has solutions of multiplicity greater than one. □ 

Remark 1. Proposition l(v) states that in general it is very unlikely for the equation i/j(z) = q to have 
multiple solutions, unless q = and E[Xi] = 0, in which case we have a double root at zero. In partic- 
ular, when doing numerical computations, we can safely assume that all solutions to ip(z) = q are simple. 

Our first goal is to identify the Mellin transform of I q , that is Ef/*" 1 ]. Our main tool is the following 
result, which we have borrowed from [10]. Note that this result is not restricted to processes with jumps 
of rational transforms. 

Proposition 2. (Verification result) Assume that Cramer's condition is satisfied: there exists z > 
such that the Laplace exponent ip(z) is finite for all z G (0, z ) and ip{9) = q for some 9 G (0, z ). If f(s) 
satisfies the following three properties 

(i) f(s) is analytic and zero-free in the strip Re(s) G (0, 1 + 9), 

(ii) f(l) = 1 and f(s + 1) = sf(s)/(q - i[>(s)) for all s G (0, 9), 
(m) I/O)!" 1 = o(exp(27r|Im(s)|)) as Im(s) -»■ oo, Re(s) G (0, 1 + 9), 
then Ef/^ 1 ] = f(s) for Re(s) G (0,1 + 9). 

Proof. The proof can be found in [10], however for sake of completeness we will present the main steps 
of the proof here. Let us denote the Mellin transform of I q as M.y{s) = Ef/*" 1 ]. First of all, the 
Cramer's condition and Lemma 2 in [20] imply that .My(s) can be extended to an analytic function 
in the vertical strip Re(s) G (0,1 + 9). Since |A^y(s)| < .My (Re(s)) we see that A^y(s) is bounded 
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in the strip Re(s) G [0/2,1 + 0/2]. It is well-known (see Lemma 2.1 in [13] or Proposition 3.1 in [4]) 
that M.y{ s ) satisfies the same functional equation as f(s). Therefore the ratio F(s) = Aiy{s)/ f(s) is 
a periodic function: F(s + 1) = F(s); and due to condition (i) F(s) can be extended to an analytic 
function in the entire complex plane. Finally, condition (iii) and boundedness of A4y{s) imply that 
F(s) = o(exp(27r|Im(s)|) in the entire complex plane, and any function which is analytic, periodic with 
period equal to one, and which satisfies this upper bound must be identically equal to a constant. Since 
F(l) = 1 we conclude that F(s) = 1, that is M Y {s) = f(s). □ 



Proposition 2 is a convenient tool which allows us to explicitly identify the Mellin transform of I q as 
a solution to the functional equation f(s + l) = sf(s)/(q — if)(s)). What makes this equation analytically 
tractable is the fact that in our case s/(q — ip(s)) is a rational function. In fact, as was pointed out by the 
referee, functional equations of the form f(s + 1) = R(s)f(s) with a rational function R(s) were studied 
by Mellin in a 1910 paper [14], where he has also investigated their connections with Gamma functions, 
Mellin transform and generalized hypergeometric functions. 

Theorem 1. Assume that q > or q = and E[Xl] < 0. For Re(s) G (0, 1 + £i) we have 

E [I- 1 ] = A 1 - x F(s) x |M, (5) 



where 



nr(i + G-«) nr(ft + ^ 

Q{s) = -p x 2=1 , (6) 



nr(i + ft-^ nr(6 + a) 

3=1 3=1 



and the constant A is defined as 



A=\, if a > 0, 



A = \fi\, if cr = and fi ^ 0, (7) 
A = q + A, if a = /i = 0. 

Proof. Let us prove this Theorem in the case a = fi = 0, the proof in other cases is similar. First of all, 
we check that the Cramer's condition in Proposition 2 is satisfied with = Let /(s) be the function 
in the right-hand side of (5). Since T(s) has no zeros in C and simple poles at s G {0, —1, —2, . . . }, it is 
clear that f(s) is analytic and zero-free in the vertical strip Re(s) G (0, 1 + 
Next, we factorize the rational function s/(q — as follows 

Hi- /;,)'"- li> • i>,Y" : 

S 3 " *~ • (8) 



g-V(s) q + X k k 

n(s-Ci) n(^+Ci) 

3=1 j=l 

One can check that this is the correct factorization, since the rational functions in the right/left hand 
sides of this equation have identical roots/poles and the same asymptotic behavior as s — > oo (recall that 
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ip(s) — > —A when s — > oo in the case a — fi — 0). Using (8) and the functional identity Y{s + 1) = sT(s) 
we check that f(s) satisfies the functional equation f(s + 1) = sf(s)/(q — ip(s)) for Re(s) G (0, 1 + CO- 
Finally, we use the asymptotic expansion 

\F(x + iy)\ = exp ^-^M + ^ - ln|j/| + 0(1)J , y ^ oo (9) 

(see formula 8.328.1 in [9]), the fact that 

j j 
M = s ^m j = K, M = s ^m j = K 

3=1 3=1 

(see Proposition l(iv)) and conclude 

= exp (-||Im(a)| + 0(ln , Im(s) oo 

in the vertical strip Re(s) G (0, 1 + Ci)- Thus we see that f(s) satisfies all three conditions of Proposition 
2, therefore (5) is true for Re(s) G (0, 1 + Ci)- □ 



Remark 2. Using Proposition 1 we see that in all cases, except when a = and fi < 0, the right- 
hand side of (5) has more gamma functions in the numerator than in the denominator. This fact and 
the asymptotic formula (9) imply that Ef/*" 1 ] decreases to zero exponentially fast as Im(s) — > oo. In the 
case a = and fi < we have the same number of gamma functions in the numerator and denominator, 
therefore (with the help of (9)) we conclude that the rate of decrease is polynomial. This implies (via 
the inverse Mellin transform) that p(x) is a smooth function on IR + , except in the case when a — and 
/i < 0. This result is consistent with Proposition 2.1 in [4] (see also Proposition 3 below). 

Formula (5) gives us the Mellin transform of I q , which uniquely characterizes the distribution of I q . 
It turns out that with little extra work we can obtain an explicit formula for the joint transform of 

Corollary 1. Assume that q > 0. For Re(tt) G (— Ci, CO an< ^ Re(s) G (0, 1 + Ci — Re(w)) we have 

E [^-'1 = ^£5 x r W x §±g. (io) 

Proof. Assume that u G (— Ci, CO anc ^ define a new measure P as the Escher transform of P 



dP 
dP 



e «X t -ii/)(u) 



Define I it) = J ' exp(X v )dv. Then we have 

oo oo 

E [e^Uf 1 ] = qj e~ qt E [e uXt /(t) 5 " 1 ] dt = q J e^'+^E [/(t) 5-1 ] dt = |e [if 1 ] , (11) 



6 



where we have denoted q = q — ip{u). Note that ip{u) < for u G (— Ci, Ci)> therefore q > 0. The Laplace 
exponent of X under the measure P is equal to ip(z) = ip(z + it) — ip(u), therefore the zeros and poles 
of the rational function ip(z) — q can be obtained from the corresponding zeros and poles of ip(z) — q by 
shifting by u units: 

Cj Cj ~ u, pj ' y pj - u, Q Cj + "• P j >-> Pj + "• 
We apply results of Theorem 1 and conclude that 

E [/?-'] = A'- x r(.) x |1±|. (12) 

Identity (10) for real values of u follows from (11) and (12), the general case follows by analytic contin- 
uation. □ 



As we have mentioned above, formula (5) uniquely characterizes the distribution of I q via the inverse 
Mellin transform 

p{x) = ^[ J E K -1 K sds - ( 13 ) 

1+iK 

In fact, p(x) can be computed explicitly, and this can be achieved in a variety of ways. One approach 
(which is quite general) is to use the fact that the right-hand side of (5) is a meromorphic function 
in C, then compute its residues, and then by shifting the contour of integration we will obtain either 
convergent series representations or complete asymptotic expansions for p(x). This is the method that 
was used in [10] in the case of hypergeometric Levy processes. The second approach, which we will follow 
in this paper, is to use the theory of the Meijer G-function. 

First of all, we want to ensure that the Mellin transform of I q given in (5) does not have multiple 
poles. One can still obtain explicit formula in the case when we do have multiple poles, however in this 
case we will have logarithmic terms and more lengthy expressions. We would like to avoid this situation, 
therefore from now on we will work under the following assumption, which ensures that T(s)Q(s) does 
not have multiple poles. 

Assumption A. 

A.l ip(z) — q has no multiple poles in the half-plane Re(z) < 0. 

A. 2 For 1 < i < M we have p t £ N. 

A. 3 For 1 < % < j < M we have pj - p { £ M. 

A. 4 i/i(z) — q has no multiple zeros in the half-plane Ke(z) > 0. 

A. 5 For l<i<j<Kwe have Q — d ^ 

Note that Assumptions A.l, A. 2 and A. 3 can be easily verified, since usually we know explicitly the 
Levy measure, and therefore the values of pj. Assumption A.l is equivalent to requiring rhj = 1 for 
j = 1, 2, . . . , J. Assumption A. 4 can be verified numerically, however due to Proposition l(v) we would 
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expect that this assumption will be satisfied in numerical experiments, unless we have a specifically 
engineered counter-example. Similarly, condition A. 5 should be satisfied in numerical experiments, since 
a "generic" polynomial equation normally does not have solutions which differ by an integer number. 
Thus Assumptions A. 4 and A. 5 are not very restrictive. 

Next, let us introduce some notations. For any w6C and n G N we define the vector [w] n G C™ as 

[w] n = [w,w, ...,w]. 

We define vectors a G C p+1 and b G as 

a = [!)[!- Plkn [1 - Pl\mi, ■■■,[!- Pj]mj, [1 + Pl] mi , [1 + P2]m 2 , + Pj]mj], 

b = [i + Ci,i + C2,...,i + a,i-Ci,i-C 2 ,---,i-C^]- 

For any x G C n and j G {1,2, . . . ,n} we will denote as (x)- 7 G C ra_1 the vector obtained from x by 
deleting the j-th coordinate, that is 

(x)^ [X\j ■ • • j •Ej— 1) •Ej+li ■ ■ ■ i •^n\- 

Finally, for x G C n and w G C, we define the vector w + x G C n as follows 

W + X = [w + Xi, W + X 2 , ■ ■ ■ ,w + x n ]. 
For u G C p and vGC 5 the generalized hypergeometric function is defined as a formal power series 

(Ui) n . . . (Up) n Z n 

(vi) n .. . {v q ) n n! ' 

where (w) n — w(w + 1) . . . (w + n — 1) is the Pochhammer symbol. Using the ratio test, it is easy to see 
that the domain of convergence of this series is (a) C when p < q (b) {\z\ < 1} when p = q + 1 (c) {0} 
when p > q + 2. In the latter case we will interpret p F q as an asymptotic expansion as z — > 0. 

Finally, let us define the following two functions (or formal power series in the case if they do not 
converge): 



F 
p 1 q 



u 






z 


V 





GAx) 



M+l 

n m-bj) n r(i 

K l<i<K i=l 

Q P+i 

i =i n ra+bj-fei) n r^-^-; 

i=K+l i=M+2 



(14) 



XX b * P+iFq_! 



1 + bj - a 

i + b, - {by 



■l) M - K x- 1 



K 



II r (% - fli) n r(i + h 



G 2 (x) = 



M+l l<i<M+l 



1=1 



p+1 



(15) 



i =1 n r(i + a l -a J ) n r( ai -6i 

i=M+2 i=Jf+l 



1 — aj + b 
1 - a, + (a)j 
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Our final result is the following proposition, which completely describes the distribution of the expo- 
nential functional. For the definition of the Meijer's G-function see formula 9.301 in [9]. 

Proposition 3. Assume that q > or q = and K[Xi] < 0. Then p(x) can be expressed in terms of 
Meijer's G-function as follows 



p(x) 



K,M+1 



-A X( 
0(1) p+1 < 



{Ax)- 1 



a 
b 



(16) 



Moreover, we have the following series representations and asymptotic expansions for p(x): 
(i) if a > then 



p(x) 
p(x) 



20(1 



G 2 ' 



20(1) 



x > 0, 
x -»■ H 



(ii) if a = and ji ^ then 



(i) if o = p, = then 



p(x) 
p(x) 



I Ml 



p(x) 



0(1) 
0(i; 



k) , 


a: > 


1/^1 


k) ) 


x < 


1/^1 



g + A 

0(1) 
g + A 

0W 



Gi((g + A)x), x ->■ +oo, 
G 2 {(q + X)x), x>0. 



Proof. Formula (16) follows from the expression for p(x) as the inverse Mellin transform (13), formula (5) 
and the definition of the Meijer's G-function (formula 9.301 in [9]). All the convergent series represen- 
tations presented in Proposition 3 follow immediately from 9.303 and 9.304 in [9], while the asymptotic 
expansions follow from Theorem 1 in [7]. □ 
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